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RESUMEN 

En el 2010 se inauguro en el Ecuador, una Presa Enrocada que es novedosa por la 
forma en que fue construida razon por la cual en este articulo se presenta en primer lugar las 
caracterfsticas tecnicas de esta mega estructura y presas similares que han sido construidas 
en el mundo para tener una mejor idea del comportamiento de las mismas. 

Luego se realiza el analisis sismico por el Metodo de Superposition Modal, para el 
sismo OBE (Operating Basis Earthquake) de un Proyecto Hidroelectrico que esta cerca de 
Mazar, utilizando dos tipos de elementos finitos, que son: el Tipo Membrana y el Q4 suavizado 
por flexion. Los resultados son comparados en terminos de desplazamientos y de fuerzas. 
Finalmente se comparan las respuestas maximas probables empleando dos criterios de 
combinacion modal. 

Palabras Claves.- Presa Enrocada. Sismo OBE. Presa del Proyecto Hidroelectrico Mazar. 
Elemento Finito Tipo Membrana. Elemento Finito Q4 suavizado por flexion. 

ABSTRACT 

In 2010 in Ecuador, a Rockfill Dam was inaugurated; it is something newfangled due to 
the manner in which it was built, for that reason in this article the technical features are first 
presented for this mega structure and similar dams that have been built worldwide in order to 
have a better idea about their behavior. 

Then the seismic analysis is carried out by means of the Modal Superposition Method, 
for the OBE (Operating Basis Earthquake) of a Hydroelectric Project that is near Mazar, using 
two types of finite elements, which are: Shell Type and Q4 element softened by flexion. The 
results are compared in terms of displacements and forces. Finally, the maximum probable 
responses are compared using two modal combination criteria. 

Keywords.- Rockfill Dam. OBE Earthquake, Dam of Mazar Hydroelectric Project, Finite 
Element of Shell Type. Q4 Finite Element softened by flexion. 


1. INTRODUCCION 

La Presa Mazar (que se abreviara PM), se halla ubicada en la region sur oriental del 
Ecuador, aproximadamente a una hora de la ciudad de Azogues, a 105 km., de la via Cuenca- 
Paute-Guarumales-Mendez, entre las cotas altitudinales de 2008 a 2400 msnm. Se encuentra 
sobre el rio Paute y fue proyectada construirse tan pronto terminara la construction de la Presa 
Paute (que se abreviara PP) para que detenga los sedimentos del rio y permita trabajar a la PP 
(tambien conocida como Presa Daniel Palacios) en optimas condiciones. 

La PM genera 160 megavatios en cambio que la PP genera 20000 megavatios. Esto se 
debe a que la primera Presa tiene como objetivo fundamental retener los sedimentos del rio 
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Paute y regular el caudal para las Presas que se hallan aguas abajo, entre ellas la PP. En la 
figura 1 se presentan datos generales de la Presa entre los que se destaca que la PM actuara 
como un reservorio para almacenar 410 millones de m 3 de agua y la PP 120 millones de m 3 
de agua; pero de esta ultima cantidad, por no construir a tiempo Mazar ya se ha almacenado 
50 millones de m 3 de sedimento. El reservorio previsto para Mazar tendra mas de 30 km., 
de largo y podra alojar sedimentos para una vida util de la Presa de 50 anos. 
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Figura 1 Ubicacion de la Pesa Mazar y datos generales. (Imagenes de Presa Mazar. Internet) 



Figura 2 Vista de Pantalla y enrocado de Presa Mazar. (Imagenes de Presa Mazar. Internet). 


Mazar es una Presa enrocada con una pantalla de hormigon, la cota de coronacion, con 
parapeto es de 2166 m., y tiene una longitud de coronacion de 340 m. A la izquierda de la 
figura 2 se presenta la Pantalla de la Presa y el agua almacenada. No se aprecia en forma 
clara el vertedero que esta a la derecha de la fotograffa izquierda. La altura maxima de la Presa 
es de 166 m., desde la cimentacion y el volumen de hormigon utilizado es de 36130 m 3 - 

A la derecha de la figura 2 se ve la Presa desde aguas abajo, se aprecia el enrocado 
con un camino de acceso a la coronacion. A la izquierda de esta fotograffa se ve el vertedero 
que sirve para regular el volumen de agua del reservorio. El nivel maximo de operacion del 
embalse es la cota 2153 m., pero en el evento de maxima crecida probable puede alcanzar la 
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cota 2165.50 m. El nivel de operation del embalse es 2098 msnm correspondiente a un 
volumnen de embalse de 101 millones de m 3 HidroPaute (2010). 

Antes de describir el vertedero se debe manifestar que el volumen del material de 
enrocado es de 5483000 m 3 . Si se compara esta cantidad con el volumen de hormigon utilizado 
en la pantalla de 36180 m 3 Se concluye que en esta Presa menos del 1% del material es 
Hormigon de tal forma que represento un gran ahorro economico, ya que el 99% es 
material de relleno. Pero tambien es importante estudiar mas sobre el probable 
comportamiento sismico que puede tener la Presa, ya que Embalses similares han presentado 
rajaduras en la pantalla de hormigon armado durante sus primeros anos de vida, sin la 
ocurrencia de ningun sismo. (INGETEC, 2007). 

En la figura 3 se muestra el vertedero de la PM durante su construction (Foto tomada 
en Agosto de 2009). La altura de cafda es de 153 m., de tal manera que entre la cota superior 
del vertedero y la cota de coronation existe una diferencia que de 13 m. Este vertedero es de 
tipo rectangular a cielo abierto sin compuertas, en el cual se utilizaron 55870 m 3 de hormigon 
armado, cantidad mayor a la empleada en la Pantalla. 



Figura 3 Vertedero de la Presa. 


En la fotografia izquierda de la figura 3 se aprecia la roca utilizada como material de 
relleno y en la fotografia de la derecha se presenta un muro que se encuentra entre la Presa, 
propiamente dicha y el vertedero. 

2. DESCRIPCION DE LA PRESA 

La pendiente de la Presa aguas arriba es 1 .405 H: 1 .0 V y aguas abajo es 1 .5 H: 1 .0 V. 
El Area A , de la cara frontal es de 48172 m 2 y la altura del embalse H es de 158 m., (sin 

considerar la cimentacion). Luego la relation A/ H 2 =1.93 este es un parametro que se 
utiliza para comparar el comportamiento de las Presas enrocadas. El espesor e de la pantalla 
de la PM se obtuvo con la siguiente ecuacion: 

e = 0.3 + 0.003// (1) 

Para la losa de borde el espesor e se hallo con la ecuacion (2). En la figura 4 se 
presentan los elementos de Hormigon Armado (HA) de la cara de la P.M. En esta figura se 
observa que la losa de borde se halla a la derecha de la pantalla junto al y al vertedero que no 
esta en la figura 4. El refuerzo de acero utilizado fue del 0.5% en direction horizontal y 
vertical. 


e = 0.6 + 0.006 H 


( 2 ) 
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En la figura 5 se presenta una seccion tfpica de la PM en el sentido transversal, en la 
que se aprecia la pantalla de HA., la cual tiene espesores de 0.77 m., en la base y de 0.30 m., 
en la corona (se encuentran con la ecuacion 1). El ancho total de la corona es 10.20 m. 
Ademas se detalla que el enrocado esta compuesto por dos tipos de materiales identificados 
como 1 y 2 en la figura 5. 

Ahora bien para el analisis sfsmico y para el diseno en general de la Presa se necesita 
conocer cual es el modulo de elasticidad E , el modulo de Poisson v y el peso especffico y . 
Para la pantalla de HA., no hay problema ya que depende de la calidad del hormigon. Sin 
embargo para el material de relleno si lo hay; inicialmente se puede pensar en extrapolar los 
datos con otras Presas que tengan similares materiales y sobre todo con una topograffa del 
canon parecida. Estos valores referenciales deberan ser ratificados o rectificados con 
simulacion numerica y ensayos de laboratorio y es asf como se lo ha hecho INGETEC (2007) 
en la PM. 

3. EXPERIENCES EN PRESAS SIMILARES 

El parametro A/ H 2 tiene como objetivo considerar la forma del encanonado de la 
Presa, en el comportamiento de un embalse enrocado, en lugar de la relacion longitud de la 
corona con respecto a la altura, LI H , debido a que se toma en cuenta el area de la cara de la 
pantalla A que esta en funcion de la topograffa del encanonado. 

INGETEC (2007) hizo un acopio de varias presas enrocadas construidas en el mundo, 
con diferentes tipos de materiales de relleno (enrocado) y vale la pena analizar por lo menos 
cual fue el asentamiento maxi mo encontrado en el proceso constructive y posteriormente con el 
paso del tiempo; cual fue el desplazamiento transversal hallado en la cara de la pantalla y el 
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modulo del material del relleno; Todo esto para entender mas sobre el comportamiento de las 
Presas Enrocadas. 


I 


28,00 


490,40 

Figura 5 Section de la Presa en direction paralela al rfo Paute. 



En la figura 6 se muestran los asentamientos (consolidation) maximos de varias Presas 
para diferentes relaciones de A/ H 2 . Aparentemente los asentamientos no dependen de esta 

relation ya que para valores de A/ H 2 bajos o altos se tienen asentamientos pequenos o 
altos. Lo que si depende es del material de relleno, ya que si se emplean materiales 
compresibles se tendran mayores asentamientos. 
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Figura 6 Asentamientos de material de relleno en Presas Enrocadas. INGETEC (2007) 


Estan subrayadas las Presas que han presentado rajaduras en la pantalla y todas 
ellas han tenido asentamientos muy altos superiores a los dos metros. Es preocupante lo 
ocurrido en la Presa Tiangshenguido en China que luego de 3 anos de operacion empezo a 
presentar rajaduras verticales en la pantalla, el asentamiento maxi mo ya esta en los 3.75 
metros. 


126 


Aguiar R., Medina A. 


Es importante analizar los datos de la Presa Campos Nova, construida en Brasil (zona 
de bajisima peligrosidad sismica) y que inicialmente presento una pequena rajadura vertical en 
la parte superior que causo cierta inquietud por las filtraciones de agua y nada mas. Pero 
cuando se tuvieron problemas en los Tuneles de Desviacion, se vieron obligados al 
desembalse completo de la Presa (Junio de 2006) y se encontro con el gran dano que muestra 
la figura 7. 





usT.: 

Figura 7 Dano en la pantalla de la Presa Campos Nova en Brasil. INGETEC (2007) 


Tabla 1 Caracteristicas generales de Presas Cam 

pos Nova y Mazar. 

Parametro 

Campos 

Nova 

Mazar 

Comentario 

Altura (m) H 

202 m. 

158 m. 

Presa de Brasil de mayor altura del 
mundo en sistema enrocado 

Pendiente aguas arriba 

1.3 H : 1.0 V 

1.4 H : 1.0 V 

Brasil no tiene problemas 
sismicos, Ecuador si. 

Pendiente aguar abajo 

1.4 H : 1.0 V 

1.5 H : 1.0 V 

Espesor en primeros 100 m. 

e=0.3+0.002H 

e=0.3+0.003H 

Mayor espesor de pantalla en 
Ecuador en los 100 m. pero luego 
el mayor espesor es de Presa de 
Brasil 

Espesor sobre los 100 m. 

e=0.3+0.005H 

e=0.3+0.003H 

Refuerzo en el centro 

0.3% en H. 
0.4 % en V. 

0.5 % en H. 
0.5 % en V. 

Mayor cuantia de acero en la 
Presa de Ecuador, en centro de 
pantalla. 

Refuerzo en el perimetro 

0.5% en H. 
0.5% en V. 

0.5% en H. 
0.5% en V. 

Igual cuantia de refuerzo. 

Material llenante en las juntas 

No 

Si 

Datos asumidos para Presa de 
Ecuador 

Refuerzo Anti-spalling 

No 

Si 

Resistencia del Hormigon 

21 MPa 

21 MPa 


Modulo de enrocado 

50-60 MPa 

150 MPa 

Con este valor se realiza el 
analisis sismico. 

Relacion Al H 2 

2.57 

1.93 


Maxima deflexion en la cara 

55 cm. 

45 cm. 

Ver figura 8. 
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En la tabla 1 se presenta las caracteristicas generales de la Presa Campos Nova de 
Brasil y Mazar. En el caso de Ecuador se han tornado mayores precauciones, como por 
ejemplo haber limitado la altura y una mejor calidad del material del enrocado, que se describe 
con mas detalle en el siguiente apartado. Sin embargo no se debe perder de vista que la 
Peligrosidad Sismica del Ecuador es una de las mas altas del mundo, a diferencia de Brasil en 
que su sismicidad es muy baja. En Aguiar (2011) se hace ver que el mega sismo de Chile de 
2010 tuvo igual magnitud que el sismo de 1906 que afecto al Sur de Colombia y Norte de 
Ecuador pero fue sentido desde Quito hasta Bogota. 

En la figura 8 se presenta en el eje horizontal, el factor de forma Al H 2 y en el eje 
vertical la relacion DIH . Siendo D la deflexion normal en la Cara y H es la altura del 
embalse. Las Presas que se hallan sobre la linea continua son aquellas que han 
experimentado dano en la Cara de la Pantalla de HA. 



Figura 8 Deflexion normal en la Cara de la Pantalla en funcion de A / H 2 INGETEC (2007) 


Para la PM el Proyectista estructural, ha identificado con linea entrecortada el valor de 
DIH que es igual a 2.85. Luego al sustituir i/ = 158m., se encuentra que la deflexion 
lateral maxima esperada en la Presa Mazar es 45 cm. 

Finalmente con relacion al Modulo Vertical del Material de relleno, se debe indicar que 
la mayor parte de Presas que presentaron dano en la Pantalla de HA., tuvieron modulos de 
elasticidad bajos como se observa en la figura 9 (se a subrayado las Presas con dano). 

Nuevamente el modulo independiente de la relacion Al H 2 , depende de la resistencia del 
material. 

Se destaca que el modulo de elasticidad del material de relleno encontrado en otras 
Presas por Plint y Filho (1998) es mayor al modulo de elasticidad vertical. Con un valor medio 
que esta entre 3 y 4. Para la Presa Enrocada de Golillas de Colombia esta relacion vale 1 .5. 
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4. PROPIEDADES DE LOS MATERIALES 

Interesa conocer cual es el modulo de Elasticidad y de Poisson de la pantalla y de los 
materiales utilizados en el enrocado en la PM, para realizar el analisis sismico. Pero antes de 
ello es importante mencionar que INGETEC calibro los modelos numericos de calculo y 
determino las Propiedades de los materiales de la Presa PM, mediante el estudio del 
comportamiento de cinco presas, que consideraron similares las mismas que se indican en la 
tabla 2 y de casos teoricos de otras Presas. 
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Figura 9 Modulo vertical del relleno en funcion de A/ H 2 INGETEC (2007) 
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Tabla 2 Caracteristicas generales de Presas con las cuales INGETEC calibro los modelos 


Presa 

Relleno 

H 

(m.) 

A/H 2 

Observaciones 

Antamina 

(Peru) 

Enrocado 

120 

5.65 

Cambio de uso de la Presa de 
embalse a retencion de material de 
mineria. Falla Bordillo en diciembre 
de 2004 (No se observa armadura 
horizontal en el dano de la Presa) 

Golillas 

(Colombia) 

Gravas 

125 

1.00 

Estribos de la Presa con inclinaciones 
superiores a 60° similar a estribo 
derecho de Presa Mazar. 

Cethana 

(Australia) 


110 

2.05 

Utilizaron e = 0.30 + 0.002 H 
Construyeron primero todo el relleno y 
posteriormente la Pantalla. (Marengo, 
2010) 

Campos Novo 
(Brasil) 

Basalto altamente 
compresible 

202 

2.57 

En tabla 1 se describe con bastante 
detenimiento la Presa. 

Cajon 

(Mexico) 

Modulo de relleno 
de 145 MPa. 

190 

3.10 

Comportamiento excelente observado 
durante funcionamiento de la Presa 

con filtraciones inferiores a 140 I/s. 


Analisis Sismico de la Presa Mazar. 


129 


A mas de ser similares es importante mencionar que se disponfa de abundante 
informacion de las Presas de la Tabla 2, en cada una de las fases de su construccion y 
despues de terminada; gracias al monitoreo de las mismas. Vale la pena estudiarlas para 
entender mas el comportamiento de este tipo de presas. 

A mas de ello INGETEC realizo ensayos de laboratorio de los materiales de la zona que 
iban a ser utilizados como material de relleno. Con toda esta informacion mas los estudios de 
las cinco Presas indicadas en la tabla 2 y de los casos teoricos; realizaron simulaciones para 
encontrar los esfuerzos y deformaciones en la pantalla de HA, las mismas que se presentan en 
la tabla 3. (Ver figura 5 para Enrocados Tipo 1 y 2) 


Tabla 3 Resultados para distintos materiales de relleno. Construccion mas llenado. INGETEC (2007) 


Modulo 

Enrocado Tipo 1 
(MPa) 

Modulo 

Enrocado Tipo 2 
(MPa) 

* 

OV 

( MPa) 

< 

(1 MPa) 

8 

(cm) 

40 

40 

20.2 

26.4 

70.5 

70 

40 

16.8 

22.0 

44.7 

70 

70 

16.0 

21.5 

44.3 

90 

40 

15.8 

20.4 

37.0 

90 

70 

14.5 

20.1 

36.2 

90 

90 

14.4 

19.5 

36.6 

110 

40 

13.9 

19.4 

31.3 

110 

70 

11.0 

18.9 

31.3 

110 

90 

14.6 

18.3 

31.3 

110 

110 

13.9 

18.1 

31.7 

150 

150 

9.9 

18.2 

25.1 

* Esfuerzo maximo horizontal 

+ Esfuerzo maximo en la direccion del talud 

8 Desplazamiento maximo en cm., perpendicular al eje de la cara. 




En base a los resultados que se han presentado en este apartado y a informacion 
encontrada en INGETEC (2007) se considero para el analisis sismico, los valores que se 
indican en la Tabla 4, en que se presenta el Modulo de Elasticidad E , de Poisson v , Peso 
especffico y . 


Tabla 4 Propiedades de los materiales con los cuales se realiza el analisis sismico. 


Material 

E 

V 

7 

(Tim 3 ) 

(MPa) 


Pantalla de HA 

20000 

0.3 

2.4 

Enrocado Tipo 1 

150 

0.3 

2.1 

Enrocado Tipo 2 

150 

0.3 

2.1 

Roca de cimentacion 

50000 

0.2 

2.4 


5. SISMO DE ANALISIS 

En el diseno sismico de Presas se consideran dos sismos denominados OBE 
(Operating Basis Earthquake) que tiene un perfodo de retorno de 140 anos y el sismo MCE 
(Maximum Credible Earthquake) con un perfodo de retorno que esta alrededor de los 970 anos. 
USCOLD (1999). Los autores de este artfculo no disponen de los espectros de diseno para la 
PM por lo que se decide trabajar con el espectro de diseno elastico OBE obtenido para una 
Presa que se encuentra cerca de Mazar. Aguiar (2010,1). Las aceleraciones espectrales estan 
definidas por las siguientes ecuaciones: 
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A e = 248 + 3720 T 
A e = 620 

A e =740.96- 366.56 T 
A e = 216.78 


T < 0.10 
0.10 <T< 0.33 
0.33 <T <1.43 
r>1.43 


(3) 


Donde A e es la aceleracion espectral elastica en gals; T es el periodo de la 

estructura. Se destaca que de acuerdo al modelo que se utilice en el analisis sismico se tendra 
“n” grados de libertad y por consiguiente “n” perfodos. Para cada periodo se hallara la 

aceleracion espectral A e . En la figura 10 se muestra el espectro de analisis. 



Figura 10 Espectro considerado en este articulo. Para el sismo OBE. 


Para el sismo OBE no se espera ningun dano en la estructura, razon por la cual se 
realiza un analisis sismico elastico considerando el Metodo de Superposicion Modal. No esta 
por demas indicar que el factor de reduccion de las fuerzas sismicas R , para este caso vale 1 . 



Figura 1 1 Grados de libertad de elemento finito considerado. 
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6. ELEMENTOS FINITOS 

Se va a realizar el analisis sismico, con un elemento finito rectangular Q4, considerando 
dos grados de libertad por nodo pero empleando dos formulaciones; uno derivado a partir del 
efecto de Membrana y el otro con el elemento Q4 suavizado por flexion. Onate (1992). En la 

figura 11 se presentan los grados de libertad del elemento finito; donde u t es la componente 
de desplazamiento en sentido horizontal y v ; es la componente de desplazamiento en sentido 
vertical; del nodo i . 

La formulacion del efecto de Membrana fue desarrollada por Guendelman (1978) y esta 
recogida en Campbell et al. (2010); la matriz (triangular superior) para las coordenadas del 
elemento indicado en la figura 1 1 , es la siguiente. 
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Donde b, a son las dimensiones del elemento finito en sentido X,Y; E es el modulo de 
elasticidad del material; v es el modulo de Poisson; t es el espesor del elemento. A los 
resultados que se obtengan con esta matriz se denominara: Metodo 1 . 

Para el Metodo 2 se trabaja con el elemento finito Q4 suavizado por flexion, cuyas 
funciones de forma son las indicadas en la tabla 5 y en la figura 12 se presentan los grados de 
libertad con los cuales se obtiene la matriz de rigidez de 12 X 12 para luego condensar a los 8 
grados de libertad indicados en la figura 1 1 . 
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Figura 12 Grados de libertad de elemento finito. 


Tabla 5 Funciones de forma de elemento finito cuadrilatero 


Coordenada 

natural 

Coordenada 

centroidal 

Deformada 

elemental 

Funcion de Forma 

1 


ul 

<fi\ =0.25(-s + lX-f + l) 

vl 

<j\ =0.25(-j + lX-f + l) 

2 

u2 

<t) 2 =0.25(s + l)(-f + l) 

v2 

+ 

1 

X 

H 

+ 

>5 
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u3 

</> 3 = 0.25(5 +ix*+i) 

v3 

(f> 3 = 0 . 25(5 +iX*+i) 

4 

u4 

<f> 4 =0.25(- s + l\t + l) 

v4 

<j> 4 =0.25(-.y + lXf + l) 


5 

u 5 

05 =l-« 2 

v5 

05 = l ~s 2 

6 

u6 

06 =1 ~t 2 

v6 

06 = 1-t 2 


En Aguiar (2010,2) se presenta en forma muy detallada el calculo de la matriz de rigidez 
del elemento finito Q4, desde la deduccion de la ecuacion (6) 

K e =\B‘ C T BdV ( 6 ) 

v 

Donde K e es la matriz de rigidez del elemento finito; B es la matriz de compatibilidad 

de deformaciones que se obtiene con las funciones de forma; C T es la matriz de elasticidad; 

dV es el diferencial de volumen. Al efectuar la integration en coordenadas naturales 
(normalizadas) la ecuacion (6) se transforma en: 

i i 

K e = ^ B‘ C T B edA 
-1 -1 


( 7 ) 
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Donde e es el espesor del elemento que se considera constante y dA es el diferencial 
de area. Para la integracion, es suficiente considerar 4 puntos de la cuadratura de Gauss. Con 
lo que se tiene: 

K e = c r B t det(J) eP, (8) 

i = 1 

Donde det(j ) es el determinante del Jacobiano que corresponde al area cooperante del punto 
que se esta considerando en la Cuadratura de Gauss; P t es el peso que esta asociado a los 
puntos de la Cuadratura de Gauss. Para cuando se consideran dos puntos de integracion (en 
cada direction) el valor de P t = 1 . 

El analisis sismico se va a realizar aplicando el Metodo de Superposition Modal pero 
utilizando dos criterios de combinacion modal, el primero es el de la Normativa de Peru (2003) 
para edificios que combina el 25% del Criterio del Valor Absoluto mas el 75% del Criterio del 
Valor Maximo Probable y el segundo criterio es cuando se emplea unicamente el Criterio del 
Valor Maximo Probable. En Aguiar (2008) estan detallados estos criterios y otros mas. 

7. ANALISIS SISMICO EN SENTIDO PARALELO AL RIO 

Se va a realizar un analisis sismico en sentido paralelo al no Paute, para ello se 
considera la seccion transversal que se indica en la figura 5 y un espesor de 1 m. La pantalla 
es de seccion variable con un ancho en la base de 0.77 m., y en la corona es de 0.30 m. 

En la figura 13 se presenta la malla de elementos finitos, se trabaja con 50 divisiones 
en sentido X y 21 divisiones en sentido Y, con lo que se tienen 1050 elementos finitos mas 21 
elementos del talon, dando un total de 1071 elementos. El ancho de la corona es de 10.20 m; 
el ancho de la base es de 468.40 m. (en la parte superior del talon); el ancho del talon en la 
parte inferior es de 40.0 m., y en la parte superior de 62 m.; la altura del talon es 20 m.; la altura 
sobre el talon es 158 m. 



En la figura 14 se muestra la discretization de las masas, las mismas que se han 
concentrado en cada nodo y se considera que los apoyos en contacto con el suelo son fijos (no 
existe desplazamiento horizontal ni vertical). 

Se realizo un analisis sismico modal, empleando el espectro indicado en la figura 10 
que esta muy bien detallado en Aguiar (2008) y se utilizo dos criterios de combinacion modal el 
de la Norma Sfsmica del Peru de 2003 y el del Valor Maximo Probable. 
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8. RESULTADOS 


En las figuras 15 y 16 se presentan los resultados hallados con el Modelo 1 y en las 
figuras 17 y 18 los resultados encontrados con el Modelo 2. Las figuras 15 y 17 muestran los 
desplazamientos, mientras que las figuras 16 y 18 indican las fuerzas. En todos estos casos se 
ha utilizado el Criterio de Combinacion modal de la Normativa Sfsmica de Peru (2003) que es 
bastante conservador. (25% del Valor Absoluto + 75% del Valor Maximo Probable). 



Figura 14 Discretization de las masas y condiciones de apoyo. 
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Figura 15 Desplazamientos maximos probables encontrados con Modelo 1. 
Criterio de Combinacion Modal de Normativa de Peru 2003. 
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Figura 16 Fuerzas maximas probables encontradas con Modelo 1. 
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Criterio de Combinacion Modal de Normativa de Peru 2003. 

A nivel de desplazamientos los resultados encontrados con los dos modelos son 
bastante parecidos. El mayor desplazamiento encontrado con el modelo 1 es 0.79 m., y con el 
modelo 2 es 0.81 m. De tal manera que a nivel de desplazamientos maximos probables los dos 
modelos reportan resultados que pueden considerarse iguales. 

Con el modelo 1 los desplazamientos son ligeramente menores a los que se hallan con 
el modelo 2. A nivel de fuerzas maxi mas en los nodos sucede lo contrario, es decir las fuerzas 
maximas halladas con el modelo 1 son mayores a las del modelo 2. A nivel de fuerzas ya 
existe una diferencia que es apreciable y se pondra mas en evidencia cuando al final del 
apartado se indique el coeficiente sfsmico. 
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Figura 17 Desplazamientos maximos encontrado con Modelo 2. 


Criterio de Combinacion Modal de Normativa de Peru 2003. 



Criterio de Combinacion Modal de Normativa de Peru 2003. 

Es importante destacar que los mayores desplazamientos laterales se presentan en el 
enrocado del material tipo 2. 

Los desplazamientos laterales en la pantalla de HA son menores a los 25 cm. Por lo 
tanto es menor al maximo permitido indicado en la tabla 3. 

En la tabla 6 se presentan los periodos de vibracion encontrados con los dos modelos, 
para los 10 primeros modos de vibracion. Los periodos de vibracion son bastante altos debido 
a que la Presa enrocada es bastante flexible en comparacion con las Presas de Hormigon. 
Pero por otra parte al tener periodos altos es beneficioso ya que se esta trabajando en la zona 
del espectro que tiene aceleraciones bajas pero desplazamientos altos. 
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El peso total de la seccion transversal analizada por metro de ancho es 81572.02 T. 
Los coeficientes sismicos que se encontraron son: 0.778 para el modelo 1 y 0.608 para el 
modelo 2. Lo que demuestra que la diferencia es notable. A mas de ello los coeficientes 
sismicos son bastante altos debido al Criterio de Combinacion Modal que es bastante 
conservador. 


Tabla 6 Periodos de vibracion encontrados con los dos modelos. 


Modo 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Periodo 
Modelo 
1 (seg.) 

2.486 

1.558 

1.056 

0.950 

0.830 

0.783 

0.760 

0.663 

0.621 

0.570 

Periodo 
Modelo 
2 (seg.) 

2.530 

1.703 

1.165 

0.954 

0.856 

0.831 

0.762 

0.679 

0.655 

0.574 


Cuando se aplica el Criterio del Valor Maximo Probable, los desplazamientos y fuerzas 
se reducen considerablemente, con relation a los obtenidos con el Criterio de la Norma 
Sfsmica de Peru (2003), que es para edificios. Se lo utilizo en Presas ya que se queria ver que 
tanto difiere con el criterio del Valor Maximo Probable. 



Figura 19 Desplazamientos maximos probables encontrados con Modelo 1. 
Criterio de Combinacion Modal del Valor Maximo Probable. 



Criterio de Combinacion Modal del Valor Maximo Probable. 

En las figuras 19 y 20 se muestran los resultados hallador con el Modelo 1 pero 
empleando el Criterio del Valor Maximo Probable. Ahora si se comparan las figuras 15 y 19 se 
aprecia que a nivel de desplazamientos existe una diferencia que esta alrededor del 15%. 
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Al comparar las figuras 16 y 20 se nota que las fuerzas obtenidas con el Criterio del 
Valor Maximo Probable disminuyeron en un 50%, aproximadamente. 


En las figuras 21 y 22 se presentan los desplazamientos y fuerzas, maxi mas probables, 
encontradas con el Modelo 2 y empleando el Criterio de Combinacion Modal del Valor Maximo 
Probable. Nuevamente se aprecia que la reduccion de los valores es notable con relacion a los 
que se encuentran con el Criterio de la Normativa sfsmica de Peru, 2003. 
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Figura 21 Desplazamientos maximos probables encontrados con Modelo 2. 
Criterio de Combinacion Modal del Valor Maximo Probable. 



Criterio de Combinacion Modal del Valor Maximo Probable. 

En la tabla 7 se muestra el coeficiente sismico hallado con los dos modelos y con los 
dos criterios de combinacion modal. Se ve que el criterio de combinacion modal que se utilice 
en el Metodo de Superposicion Modal influye mucho mas que la comparacion entre modelos de 
elementos finitos que se utilicen. 

En base a lo obtenido, definitivamente se debe descartar el criterio del 25 % del Valor 
Absoluto mas el 75 % del Valor Maximo Probable, ya que es demasiado conservador y lo que 
se debe utilizar es simplemente el Criterio del Valor Maximo Probable. 

Con este ultimo criterio se aprecia que los resultados encontrados a nivel de fuerzas y 
desplazamientos, con los dos modelos son bastante parecidos. 


Tabla 7 Coeficientes sismicos hallados con dos modelos y dos criterios de combinacion modal. 


Criterio de Combinacion Modal 

Modelo 1 

Modelo 2 

Norma Tecnica de Peru 

0.778 

0.608 

Valor Maximo Probable 

0.311 

0.279 
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9. COMENTARIOS Y CONCLUSIONES 

Llamo mucho la atencion la construction de una Presa Enrocada por el ahorro 
considerable ahorro de dinero que se tiene en relacion a una Presa de Hormigon Rodillado de 
similares caracteristicas. Por este motivo es muy importante estudiar este nuevo sistema 
constructive para el Ecuador y ver las experiencias que han tenido Presas similares en otras 
partes del mundo. Este es un primer articulo y se aspira continuar investigando siempre y 
cuando HIDROPAUTE Y CONELEC nos proporcionen la informacion sismica que ha sido 
registrada en el Proyecto Mazar. 

Se ha realizado el analisis sismico elastico de la Presa Enrocada del Proyecto Mazar, 
para un sismo OBE, mediante el Metodo de Superposition Modal; utilizando dos tipos de 
elementos finitos y dos criterios de combinacion modal, del estudio realizado se desprenden las 
siguientes conclusiones. 

• El criterio de combinacion modal que combina el 25% del Valor Maximo Absoluto mas 
el 75% del Criterio Valor Maximo Probable, es demasiado conservador y no se lo 
puede utilizar en Presas. Este es un criterio adecuado para edificios. Para Presas lo 
mejor es utilizar el criterio del Valor Maximo Probable, con relacion al anterior. 

• A nivel de desplazamientos nodales el elemento Finito con efecto Membrana reporta 
resultados muy parecidos a los que se hallan con el elemento Finito Q4 de tal manera 
que no se debe descartar el uso del Modelo 1 en el analisis sismico de Presas, al 
contrario al efecto Membrana se debe anadir el efecto de los giros y comparar los 
resultados que se hallen con el elemento Finito Q4. 
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